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Abstract 


This  report  presents  three  different  approaches  for  determining  exceedence  probabilities 
of  ship  sea  loads  using  linear  strip  theory.  The  first  two  approaches  use  nominal  maximum  loads 
for  hourly  and  variable  duration  seaways,  while  a  third  approach  considers  individual  load  cycle 
amplitudes  in  all  seaways.  Of  the  three  approaches,  the  load  cycle  approach  appears  to  be  most 
useful  and  can  be  used  for  both  fatigue  and  ultimate  load  computations.  Sample  computations 
for  the  Canadian  Patrol  Frigate  demonstrate  the  application  of  the  methods.  Gumbel  distribu¬ 
tions  provide  good  fits  to  computed  lifetime  load  exceedence  probabilities.  For  load  cycle  and 
nominal  hourly  maximum  loads,  Weibull  distributions  provide  superior  fits.  Future  work  should 
include  three-dimensional  hydrodynamic  forces,  nonlinear  load  effects,  and  the  influence  of  wave 
conditions  on  ship  speed  and  heading. 


Resume 

Le  present  rapport  fait  etat  de  trois  approches  differentes  utilisees  pour  determiner, 
a  1’aide  d’une  theorie  de  depouillement  lineaire,  l’indice  de  depassement  des  charges  exercees 
par  la  mer  sur  un  navire.  Les  deux  premieres  approches  reposent  sur  le  calcul  des  charges 
maximales  nominales  pour  etat  de  mer  de  duree  variable  et  de  duree  d’une  heure,  alors  que  la 
troisieme  approche  se  fonde  sur  l’amplitude  de  cycles  de  charge  individuels  dans  des  mers  de 
tous  les  etats.  Parmi  les  trois  approches,  celle  des  cycles  de  charge  semble  etre  la  plus  utile  et 
peut  servir  tant  aux  calculs  de  fatigue  qu’aux  calculs  de  la  charge  de  rupture.  Des  echantillons 
de  tels  calculs  effectues  sur  une  ffegate  canadienne  de  patrouille  demontrent  l’application  de 
cette  methode.  Les  distributions  de  Gumbel  donnent  une  bonne  equivalence  des  probability  de 
depassement  des  charges  de  duree  calculee.  En  ce  qui  concerne  le  cycle  de  charge  et  les  charges 
maximales  horaires  nominales,  les  distributions  Weibull  donnent  des  equivalences  superieures. 
Les  prochains  travaux  devraient  examiner  les  forces  hydro dynamiques  en  trois  dimensions,  les 
effets  de  charges  non  lineaires  et  l’influence  de  l’etat  des  vagues  sur  la  vitesse  et  la  direction  du 
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navire. 


DREA  TM/98/215 

PREDICTION  OF  EXTREME  AND  FATIGUE  SEA  LOADS 
TJSTNG  LINEAR  THEORY 

by 

Kevin  A.  McTaggart 
EXECUTIVE  SUMMARY 


Introduction 

The  rational  design  of  ship  structures  requires  application  of  probabilistic  methods.  This 
technical  memorandum  presents  three  different  methods  for  computing  exceedence  probabilities 
for  fatigue  and  lifetime  maximum  sea  loads.  The  first  approach  determines  a  nominal  maximum 
sea  load  for  each  hourly  seaway  that  the  ship  encounters.  The  second  approach  computes  the 
nominal  maximum  load  for  each  seaway  based  on  its  expected  encounter  duration  during  the 
ship  life.  The  third  approach  considers  the  variation  of  individual  load  cycle  amplitudes  for 
all  conditions.  Sample  computations  for  the  Canadian  Patrol  Frigate  (CPF)  demonstrate  the 
application  of  the  three  methods. 

Principal  Results 

The  nominal  hourly,  nominal  seaway,  and  load  cycle  approaches  all  give  useful  results;  how¬ 
ever,  the  load  cycle  approach  appears  to  be  the  most  rigorous  and  useful.  Output  from  the  load 
cycle  approach  includes  conditional  distributions  for  ship  and  wave  conditions  when  sea  loads 
exceed  a  specified  threshold  value.  Gumbel  distributions  provide  good  fits  to  computed  lifetime 
load  exceedence  probabilities.  For  load  cycle  and  nominal  hourly  maximum  loads,  Weibull  dis¬ 
tributions  provide  superior  fits.  For  extreme  lifetime  loads  acting  on  the  CPF,  the  midships 
vertical  bending  moment  is  approximately  three  times  greater  than  the  midships  horizontal 
bending  moment. 

Significance  of  Results 

It  is  recommended  that  the  load  cycle  approach  be  used  to  determine  sea  load  exceedence 
probabilities  for  both  fatigue  and  lifetime  extreme  computations.  Conditional  probabilities  of 
ship  and  sea  conditions  given  high  sea  loads  can  be  used  to  determine  which  conditions  warrant 
more  detailed  analysis.  The  current  sample  computations  for  the  CPF  indicate  that  maxi¬ 
mum  lifetime  loads  are  likely  to  be  caused  by  significant  wave  heights  of  approximately  15  m, 
suggesting  that  nonlinear  effects  will  be  very  important. 

Future  Plans 

The  current  load  cycle  approach  will  likely  be  enhanced  in  several  ways.  A  program  will 
be  developed  to  determine  composite  wave  scattergrams  based  on  areas  and  seasons  of  ship 
operations.  Future  developments  should  consider  the  dependence  of  ship  operations  on  wave 
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conditions.  For  example,  a  ship  will  usually  reduce  speed  in  severe  head  seas  to  reduce  the  inci¬ 
dence  of  slamming.  Although  the  current  implementation  is  based  on  strip  theory,  it  should  be 
extended  to  include  three-dimensional  hydrodynamic  effects,  which  are  important  for  the  CPF. 
Nonlinear  effects  will  be  important  for  severe  sea  conditions.  DREA  is  currently  sponsoring 
development  of  time  domain  methods  which  include  nonlinear  forces.  Attention  should  also  be 
given  to  accurate  statistical  modelling  of  extreme  seas  which  influence  ultimate  strength  design. 


IV 


Contents 

Abstract  ii 

Executive  Summary  iii 

Table  of  Contents  v 

Notation  vi 

1  Introduction  1 

2  Nominal  Hourly  Load  Approach  1 

3  Nominal  Seaway  Load  Approach  3 

4  Load  Cycle  Approach  4 

5  Fitting  of  Gumbel  and  Weibull  Distribution  Parameters  5 

6  Numerical  Implementation  6 

7  Example  Calculations  for  Canadian  Patrol  Frigate  7 

7.1  Vertical  Bending  Moment  at  Midships .  8 

7.2  Horizontal  Bending  Moment  at  Midships .  14 

8  Recommendations  for  Future  Work  20 

8.1  Wave  Scattergrams  for  Ship  Missions .  20 

8.2  Correlations  between  Wave  Conditions  and  Ship  Operations .  20 

8.3  Three  Dimensional  Hydrodynamic  Effects . 20 

8.4  Nonlinear  Wave  Loads .  20 

8.5  Wave  Data  for  Large  Wave  Heights .  21 

9  Conclusions  21 

Appendices  22 

A  EXTRMLIN  Input  22 

A.l  Input  Records . 22 

A.  2  Sample  Input .  25 

B  EXTRMLIN  Output  26 

B. l  Description  of  Output .  26 

B.2  Sample  Output .  27 

References  45 


v 


Notation 


D 

Fx(X) 
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XCycle 

X* 
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Xi 
X, 


seaway, a 


a 

ax 

P 

f? 

Pi 

Px 

7 

5 

Px 

a 

ox 


duration 

cumulative  distribution  function  of  variable  X 

probability  density  function  of  variable  X 

zero-crossing  frequency 

mean  zero-crossing  frequency 

significant  wave  height 

lower  bound  of  significant  wave  height  range 

Hs, lower  value  for  load  event  with  rank  i 

middle  value  of  significant  wave  height  range 

Weibull  distribution  parameter 

number  of  load  cycles 

number  of  hours  of  ship  operation 

number  of  seaways  encountered  during  ship  life 

number  of  discretized  values  of  variable  X 

probability  of  occurrence  for  discrete  value  Xi 

exceedence  probability  of  variable  X 

conditional  exceedence  probability  of  variable  X  given  Y 

peak  wave  period 

Tp  value  for  load  event  with  rank  i 

zero-crossing  period 

Gumbel  distribution  parameter  for  variable  X 
ship  speed 

ship  speed  for  load  event  with  rank  i 
discretized  value  of  Vs 
random  variable 
load  cycle  amplitude 

minimum  threshold  for  load  cycle  amplitude 

nominal  maximum  with  exceedence  probability  a  over  duration  D 
Xn,a  value  for  load  event  with  rank  i 
discretized  value  of  random  variable  X 

nominal  maximum  with  exceedence  probability  a  in  a  seaway 
exceedence  probability  for  nominal  hourly  maximum  load 
Gumbel  distribution  parameter  for  variable  X 
incident  sea  direction  (relative  to  ship  speed) 

P  value  for  load  event  with  rank  i 
discretized  value  of  (3 

Weibull  distribution  parameter  for  variable  X 
Euler’s  constant  (~  0.5772) 
lifetime  exceedence  probability 
mean  value  of  variable  X 

standard  deviation  (RMS  value)  of  random  process 
standard  deviation  of  variable  X 
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1  Introduction 


A  major  trend  in  current  ship  structural  design  and  maintenance  is  the  evolution  toward 
probabilistic  approaches.  Mansour  et  al.  [1]  give  an  overview  of  the  application  of  reliability 
methods  to  ship  structures.  Within  DND,  the  Improved  Ship  Structural  Maintenance  Manage¬ 
ment  (ISSMM)  project  is  developing  probabilistic  approaches  for  naval  ships. 

The  hydrodynamic  loads  on  ship  structures  axe  highly  variable;  thus,  any  probabilistic  ap¬ 
proach  to  ship  structures  must  include  a  probabilistic  treatment  of  the  hydrodynamic  loads. 
A  rational  probabilistic  approach  can  determine  the  statistical  distribution  of  sea  loads,  which 
can  then  be  used  as  input  for  a  structural  reliability  model.  Several  authors  present  approaches 
for  determining  statistical  distributions  of  sea  loads.  Sikora  et  al.  [2]  present  a  semi-empirical 
approach  for  estimating  vertical  bending  moments  including  nonlinear  and  whipping  effects. 
Guedes  Soares  and  Schellin  [3]  give  a  rational  method  for  including  nonlinear  effects  when  deriv¬ 
ing  long-term  distributions  for  vertical  bending  moment.  Jensen  and  Dogliani  [4]  use  a  similar 
approach  based  on  quadratic  strip  theory  for  predicting  bending  moment  distributions.  They 
indicate  that  nonlinear  effects  can  be  significant  for  ships  that  are  not  wall-sided. 

The  current  report  describes  initial  efforts  within  DND  to  derive  statistical  distributions 
of  sea  loads  acting  on  a  naval  frigate.  The  computer  program  SHIPM07  [5]  provides  linear 
predictions  of  sea  loads.  SHIPM07  provides  generally  good  predictions  of  sea  loads  [6,  7],  with 
some  overprediction  (typically  20-30  percent)  of  vertical  bending  moment  for  a  transom  stern 
frigate  such  as  the  Canadian  Patrol  Frigate  (CPF).  The  linear  sea  load  predictions  of  SHIPM07 
are  coupled  with  wave  statistics  from  British  Maritime  Technology  (BMT)  [8].  Three  different 
methods  are  presented  for  determining  load  exceedence  probabilities  during  the  life  of  a  ship.  The 
first  method  is  based  on  the  nominal  maximum  loads  in  all  of  the  one  hour  seaways  encountered 
by  the  ship.  The  second  method  is  based  on  the  nominal  maximum  loads  using  the  expected 
duration  of  each  seaway  during  the  lifetime  of  the  ship.  The  third  method  considers  all  of  the 
individual  wave  load  cycles  encountered  by  the  ship. 

The  following  section  presents  the  theory  for  the  nominal  hourly  load  approach,  which  is 
followed  by  the  nominal  seaway  load  approach.  Section  4  gives  the  alternative  approach  based 
on  individual  load  cycles.  The  fitting  of  Weibull  and  Gumbel  distributions  to  derived  data 
is  presented  in  Section  5.  The  numerical  implementation  of  the  theory  is  given  in  Section  6. 
Section  7  gives  an  example  of  application  of  the  methods  to  the  CPF.  Recommendations  for 
future  work  to  include  nonlinear  effects  and  mission  profiles  are  given  in  Section  8,  which  is 
followed  by  final  conclusions. 


2  Nominal  Hourly  Load  Approach 

In  the  first  approach  for  computing  sea  load  probabilities,  a  nominal  hourly  maximum  load 
will  be  taken  for  each  combination  of  ship  speed,  heading,  significant  wave  height,  and  peak  wave 
period.  Assumptions  of  linearity  and  narrow  bandedness  permit  load  amplitude  in  a  random 
seaway  to  be  modelled  using  a  Rayleigh  distribution  as  follows: 


fx cycle  (X cycle)  =  eXp 


_  v'2 
^  cycle 


2a2 


(2.1) 
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where  fxcycu  (X cycle)  is  the  probability  density  function  for  load  amplitude  Xcyde  in  a  single  load 
cycle  and  a  is  RMS  load  in  a  seaway.  The  nominal  maximum  load  XjjtCl  with  probability  of 
exceedence  a  during  duration  D  is  given  by: 


(2.2) 


where  Tz  is  the  zero-crossing  period  of  the  load  process.  For  the  present  analysis,  the  duration  D 
is  taken  as  being  one  hour  because  a  loading  process  can  generally  be  considered  stationary  (i.e. 
the  standard  deviation  a  and  zero  crossing  period  Tz  will  not  change  significantly)  during  a  one 
hour  period.  Table  1  shows  the  ratio  of  maximum  to  RMS  load  for  nominal  seaway  conditions. 
For  the  example  given,  the  nominal  maximum  load  increases  by  approximately  ten  percent  as 
the  exceedence  probability  is  decreased  by  a  factor  of  ten.  For  structural  design,  an  exceedence 
probability  of  the  order  of  0.01  would  likely  be  acceptable  for  computing  the  nominal  maximum 
load  in  a  one  hour  seaway. 


Table  1:  Ratio  of  Maximum  to  RMS  Load  for  One  Hour  Duration  and  Zero-Crossing  Period  of 
10  Seconds 


Exceedence  probability  a 

X  hour, a, 

0.1 

4.05 

0.01 

4.58 

0.001 

5.06 

0.0001 

5.49 

The  concept  of  the  nominal  hourly  load  in  a  seaway  simplifies  the  computation  of  long  term 
load  exceedence  probabilities.  The  current  analysis  assumes  that  ship  speed  and  heading  are 
independent  variables  but  considers  the  joint  distribution  of  significant  wave  height  and  peak 
wave  period,  which  can  be  obtained  from  wave  climate  databases  such  as  BMT  Global  Wave 
Statistics.  The  procedure  for  computing  long  term  load  exceedence  probabilities  is  as  follows: 

1.  compute  reference  RMS  ship  loads  and  zero-crossing  periods  in  irregular  seaways  for  all 
relevant  combinations  of  ship  speed  Vs,  heading  /?,  and  peak  wave  period  Tp  (only  one 
significant  wave  height  is  required  per  peak  wave  period), 

2.  provide  input  probabilities  for  ship  speed  and  heading,  and  a  wave  statistics  file  with  joint 
distribution  of  wave  height  and  peak  wave  period, 

3.  compute  nominal  hourly  wave  loads  for  all  existing  combinations  of  ship  speed,  heading, 
peak  wave  period,  and  significant  wave  height, 

4.  rank  nominal  maximum  wave  loads  from  highest  to  lowest, 

5.  compute  long  term  exceedence  probabilities  for  nominal  maximum  wave  loads. 


2 


When  performing  step  3  above,  the  nominal  hourly  load  Xhour,a(Vs)  P,  Hs, lower,  Tp)  is  computed 
at  the  lower  wave  height  H$ti0 wer  for  each  wave  height  range  in  the  wave  scattergram.  The 
exceedence  probability  QxhoUT<a  for  each  of  the  ranked  loads  is  then  computed  as  follows: 

.!?))  = 

Qx*„„r,„  +  p(Vi) pm p(,Hiilcw„,rp)  (2.3) 

where  the  superscript  i  on  each  variable  denotes  the  value  of  the  variable  associated  with  the 
load  event  having  rank  i,  with  i  =  1  denoting  the  highest  long  term  load.  A  typical  computation 
procedure  will  have  approximately  7  ship  headings,  15  significant  wave  heights,  and  11  peak 
wave  periods,  although  many  of  the  combinations  of  wave  height  and  period  can  be  neglected 
because  they  will  not  occur. 

When  considering  representative  wave  periods  in  irregular  seas,  it  should  be  noted  that 
SHIPM07  uses  peak  wave  period  Tp,  while  BMT  Global  Wave  Statistics  uses  zero-crossing  wave 
period  Tz.  The  following  conversion  can  be  used: 

Tp  =  1.408  Tz  (2.4) 

A  key  benefit  of  the  nominal  hourly  load  approach  is  that  it  applies  rankings  to  discrete 
combinations  of  ship  speed,  heading,  significant  wave  height,  and  peak  wave  period;  thus,  the 
approach  indicates  which  conditions  produce  the  highest  loads.  This  information  can  form  the 
basis  for  more  sophisticated  treatment  of  nonlinear  effects  and  possible  correlation  between 
variables  such  as  significant  wave  height  and  ship  speed. 

The  current  procedure  estimates  exceedence  probabilities  for  maximum  hourly  loads.  The 
following  equation  can  be  used  to  predict  the  probability  distribution  of  the  nominal  maximum 
load  during  ship  life: 

Fx„,(Xui,)  =  [Fx,.,,a(A:tar,„)]''w  (2.5) 

where  Fxlife(Xufe)  is  the  cumulative  distribution  function  of  the  lifetime  maximum  load  Xufe 
and  Nhour  is  the  number  of  hours  of  operation  during  the  ship  life.  For  a  ship  spending  30 
percent  of  its  30  year  life  at  sea,  the  number  of  hours  at  sea  will  be  approximately  80,000. 

3  Nominal  Seaway  Load  Approach 

Equation  (2.2)  given  above  indicates  that  the  maximum  load  with  exceedence  probability  a  in 
a  seaway  is  dependent  upon  seaway  duration  D.  When  considering  lifetime  maximum  loads,  the 
above  approach  does  not  consider  that  different  seaways  will  have  different  expected  durations, 
and  this  may  affect  the  maximum  loads.  To  circumvent  this  problem,  nominal  maximum  loads 
can  be  considered  on  the  basis  of  seaway  rather  than  on  the  basis  of  hour. 

For  the  present  discussion,  a  seaway  is  considered  to  be  a  combination  of  ship  speed,  ship 
heading,  significant  wave  height,  and  peak  wave  period.  The  number  of  different  seaways  en¬ 
countered  by  the  ship  is  given  by: 

N seaway  =  %  x  Np  x  Number  of  cells  with  p(Hsjower,Tp)  >  0  (3.1) 
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The  average  duration  for  each  seaway  is  given  by: 

D(Vs,frHs,Tp)  =  Lp(at  sea)  p(V$)p(f3)p(Hs,Tp)  (3.2) 

where  L  is  the  ship  life  and  p(at  sea)  is  the  fraction  of  the  ship  life  spent  at  sea.  The  nominal 
maximum  load  for  each  seaway,  Xseaway >a,  is  calculated  using  Equation  (2.2)  using  the  seaway 
duration  from  Equation  (3.2).  The  exceedence  probabilities  for  nominal  seaway  maximum  loads 
are  evaluated  using  the  same  procedure  as  for  nominal  hourly  maximum  loads.  The  following 
equation  gives  the  cumulative  distribution  for  lifetime  maximum  load  based  on  the  nominal 
seaway  maximum  loads: 

FXiife{Xlife)  =  [^Xseaway (XSeaway,a)]  v  (3-3) 


4  Load  Cycle  Approach 


The  nominal  hourly  and  seaway  load  approaches  developed  above  have  two  main  shortcom¬ 
ings.  The  use  of  a  nominal  maximum  load  does  not  rigorously  consider  the  statistical  variation 
of  maximum  load  during  each  hour  or  seaway.  The  second  shortcoming  is  that  the  nominal  load 
approaches  do  not  provide  the  statistical  distribution  of  individual  load  cycle  amplitudes,  which 
is  typically  used  for  fatigue  computations. 

Exceedence  probabilities  for  load  cycle  amplitudes  can  be  computed  using  the  procedure 
described  by  Paulling  [9].  Again  assuming  a  Rayleigh  distribution  for  load  cycle  amplitude  in  a 
given  seaway,  the  long  term  exceedence  probability  for  load  cycle  amplitude  in  all  seaways  is: 


Q X cycle  (X cycle  )  — 

Nys  Nj3  Nh$  XTp 


1VVS  iV/3  "Tp  / 

£  £  £  £  piy.  —i)  P{Pj )  P{Hs,mid—k}  /)  -7  exp  I 

1=1  j=l  k= 1  1=1  JZ  \ 


_ v 2  ' 

~~cyde 

2  a2 


(4.1) 


where  X^de  is  load  amplitude  in  a  single  loading  cycle,  HS  Jfll(i_k  is  the  middle  value  of  significant 
wave  height  for  range  k,  fz  is  zero-crossing  frequency  in  a  seaway,  and  fz  is  average  long-term 
zero-crossing  loading  frequency  given  by: 


NVs  Np  NHs  NTp 

/z  =  ELE  £  P(VS^)  pifr)  p(Hs<mid-k,Tp^i)  fz  (4.2) 

i= 1  j=l  k=  1  Z=1 

The  probability  distribution  of  maximum  lifetime  load  can  then  be  computed  as  follows: 

Fxlife(Xlife )  =  [Fxcycle{Xoycle)]NcyCU  (4.3) 

where  N^de  is  the  number  of  load  cycles  in  the  ship  life,  which  is  given  by: 

Ncycle  =  — V^gUr  (4.4) 

J  z 

A  naval  frigate  will  experience  of  the  order  of  107  wave  loading  cycles  during  its  life. 

When  designing  for  extreme  loads,  it  is  useful  to  know  which  conditions  cause  those  extreme 
loads.  For  loads  exceeding  a  threshold  level  X ^e,  the  theorem  of  Bayes  can  be  applied  to 
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determine  conditional  probability  distributions  for  ship  speed,  heading,  and  wave  conditions. 
For  example,  the  following  equation  gives  the  conditional  probability  distribution  of  ship  speed 
for  load  amplitudes  exceeding  X*yde : 


P(V,- t\Xcyde>X"clde)  = 


QXcycle  (X cycle) 


(4.5) 


Similar  equations  can  be  applied  to  determine  conditional  distributions  for  ship  heading,  wave 
height,  and  wave  period. 


5  Fitting  of  Gumbel  and  Weibull  Distribution  Parameters 


Once  load  probabilities  are  computed  using  any  of  the  three  approaches  described  above,  it 
is  useful  to  determine  fitted  statistical  distributions  for  load  probabilities.  A  fitted  statistical 
distribution  allows  load  probabilities  to  be  described  by  a  small  number  of  parameters.  A  fitted 
distribution  also  facilitates  interpolation  or  extrapolation  of  load  probabilities. 

Thoft-Christensen  and  Baker  [10]  give  a  useful  overview  of  statistical  distributions  used  in 
engineering  reliability.  For  sea  loads  acting  on  a  ship,  the  Gumbel  and  Weibull  distributions  are 
commonly  used.  The  cumulative  distribution  function  of  the  Gumbel  distribution  is  as  follows: 


Fx(X)  =  exp  {—  exp  [ax(X  —  «x)]}  (5.1) 

where  ux  and  ax  are  distribution  parameters.  The  mean  fix  and  standard  deviation  ax  of  the 
Gumbel  distribution  are  related  to  the  distribution  parameters  as  follows: 


fix  =  ux  +  — 

ax 

ax  = 

where  7  is  Euler’s  constant  («  0.5772). 

The  Weibull  distribution  can  be  described  by  2  or  3  parameters,  where  the  third  parameter 
represents  a  lower  threshold  (typically  non-zero)  for  the  variable  X.  The  present  study  considers 
the  2-parameter  distribution,  for  which  exceedence  probability  is  given  by: 


(5.2) 

(5.3) 


Qx(X ) 


exp 


'x_yx 

v kx  ) 


(5.4) 


where  kx  and  0x  are  distribution  parameters.  The  mean  and  standard  deviation  for  the  2- 
parameter  Weibull  distribution  are: 


px  =  kx  r  (1  +  J- 

n  -  kxHl  +  l)-T2(l  +  ^ 

where  T(y)  denotes  the  Gamma  function  with  argument  y. 


(5.5) 

(5.6) 
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Two  methods  are  available  for  determining  fitted  Gumbel  and  Weibull  distribution  parame¬ 
ters  for  sea  loads.  The  method  of  moments  relates  the  distribution  parameters  to  the  computed 
mean  and  standard  deviation  for  the  variable  X.  The  second  method  uses  a  least  squares  linear 
fit  between  the  variable  X  or  its  transformation  (e.g.  lnX)  and  exceedence  probability  Qx(X) 
or  its  transformation  (e.g.  lnQxCA)).  A  significant  advantage  of  the  least  squares  method  is 
that  it  can  be  applied  to  a  limited  probability  range  of  greatest  interest  for  design,  such  as  the 
upper  range  of  sea  loads  or  the  lower  range  of  ship  strength.  Applying  a  least  squares  fit  to  a 
limited  range  of  interest  generally  produces  a  better  statistical  fit  within  that  range. 

The  equation  for  determining  a  least  squares  fit  of  Gumbel  parameters  is  as  follows: 

X  =  —\n[-ln(Fx(X))}  +  ux  (5.7) 

ax 

Placing  the  variable  X  on  the  left-hand  side  of  the  above  equation  minimizes  the  error  in  X  for 
a  given  value  of  In  [—  ln(jFx(.X’))].  The  regression  fit  has  a  slope  of  1  /ax  and  an  intercept  of 
ux-  For  the  Weibull  distribution,  the  parameters  can  be  determined  through  a  least-squares  fit 
of  the  following: 


InX  =  +  lnkx  (5.8) 

PX 

The  resulting  regression  fit  has  a  slope  of  1  /fix  and  an  intercept  of  In  kx-  When  evaluating 
Equations  (5.7)  and  (5.8),  the  terms  In  [— ln(Fx(A))]  and  In  [— ln(Qx(X))]  become  undefined 
when  the  exceedence  probability  Qx(X)  approaches  either  zero  or  one;  therefore,  a  Gumbel 
or  Weibull  least  squares  fit  cannot  include  points  with  exceedence  probabilities  of  zero  or  one. 
Furthermore,  the  term  In  X  in  Equation  (5.8)  becomes  undefined  when  X  approaches  zero. 
When  evaluating  lateral  sea  loads,  the  occurrence  of  zero  values  in  head  and  following  long 
crested  seas  necessitates  care  when  applying  Gumbel  and  Weibull  fits. 

6  Numerical  Implementation 

The  computer  program  EXTRMLIN  implements  the  nominal  hourly  load,  nominal  seaway 
load,  and  load  cycle  approaches  described  above.  Appendix  A.l  gives  the  input  format. 

Numerical  implementation  was  relatively  straightforward.  The  most  critical  numerical  as¬ 
pect  is  the  accurate  computation  of  probabilities.  All  probabilities  are  stored  using  double 
precision,  which  is  required  because  exceedence  probabilities  of  interest  range  between  10~10 
and  1.0.  The  input  wave  scattergram  is  normalized  to  correct  any  roundoff  errors.  When  com¬ 
puting  probabilities  for  sea  loads,  values  associated  with  small  exceedence  probabilities  are  of 
greatest  interest.  To  minimize  roundoff  errors  for  exceedence  probabilities,  most  probability 
computations  are  done  using  exceedence  probability  Qx(X )  rather  than  the  cumulative  distri¬ 
bution  function  Fx(X).  The  relationship  between  the  exceedence  probability  and  cumulative 
distribution  function  is: 


Qx(X)  =  1  -  Fx(X)  (6.1) 

Equations  2.3  and  4.1  give  examples  of  formulations  where  exceedence  probability  is  used  to 
maintain  accuracy. 
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When  determining  lifetime  maxima  from  large  numbers  of  events  using  Equations  2.5  and 
4.3,  numerical  precision  becomes  a  problem  for  small  exceedence  probabilities  because  of  inaccu¬ 
racies  in  Fx(X)  as  it  approaches  unity.  When  determining  the  cumulative  distribution  function 
Fxmax  n(Xmax,n)  of  the  maximum  of  n  samples  of  variable  X,  the  following  equation  can  help 
to  preserve  numerical  accuracy: 

In  FxmaXin  (Xmax.n)  =  n\nFx(X)  (6.2) 

For  small  values  of  Qxmax,n(X™ax,n),  the  following  Taylor  series  expansion  can  be  used  to  main¬ 
tain  numerical  accuracy  when  evaluating  exceedence  probability  from  In Fxmax<n {Xmax,n)'- 


!  -xr  IT-1  /  -y  \  [In  FXmax,n  (Xmax,n )  j 

Qxmax,n\-X-max,n)  ~  ~  ln  rXmax,n  \^-max,n)  ~  2 

[In  FXmax<n  (XmaXtn)}  [In  FxmaXin  (^mai,n)] 
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(6.3) 


For  Qxmaxn(Xmax,n)  <  0.01,  Equation  (6.3)  gives  a  relative  error  of  less  than  one  millionth 
of  the  exceedence  probability  (e.g.  less  than  10-8  for  Qxmax,n(Xmax,n)  =  0.01).  The  following 
Taylor  series  expansion  can  simplify  the  evaluation  of  In FxmaXin{Xmax,n)' 


111  F Xmax,n  {Xmax,n  )  ~ 


QxiX)  1  (Qx(X))2 
Fx(X)  +  2  \Fx(X)J 


lfQx(X))3 

3\Fx(X)J 


1  fQx(X))4 

4  \Fx(X)J 


(6.4) 


The  above  series  will  give  a  relative  error  in  In Fx^^^iXmax^n)  of  less  than  a  millionth  when 
Qx{X)  is  less  than  0.01. 


7  Example  Calculations  for  Canadian  Patrol  Frigate 

Sample  calculations  have  been  performed  for  the  Canadian  Patrol  Frigate  based  on  the 
sample  case  of  the  SHIPM07  manual  [5].  Appendix  A.2  gives  sample  input  for  the  program 
EXTRMLIN  and  Appendix  B.2  gives  sample  output.  This  section  presents  exceedence  proba¬ 
bilities  for  vertical  and  horizontal  bending  moments  at  midships. 

For  the  sample  case,  the  ship  has  a  30  year  life  and  spends  30  percent  of  its  time  at  sea.  All 
time  at  sea  is  spent  in  sea  area  15  [8],  which  extends  from  Nova  Scotia  and  Newfoundland  to 
approximately  halfway  across  the  Atlantic.  This  area  has  a  relatively  severe  wave  climate;  thus, 
the  resulting  sea  loads  are  likely  conservative.  The  ship  divides  its  time  evenly  between  travelling 
at  10  and  18  knots.  Waves  are  assumed  to  be  long  crested  and  modelled  by  a  Bretschneider 
spectrum,  with  a  uniform  distribution  of  relative  sea  directions. 

For  fitting  of  Weibull  and  Gumbel  distributions,  the  input  range  of  exceedence  probabilities 
is  0.0  <  Qx{X)  <  0.9.  Trial  computations  with  the  full  range  of  exceedence  probabilities 
0.0  <  Qx(X)  <1.0  gave  inferior  matching  between  computed  and  fitted  distributions.  Initial 
fitting  of  distribution  parameters  to  computed  exceedence  probabilities  for  vertical  bending 
moment  was  relatively  simple.  Additional  fitting  of  parameters  for  horizontal  bending  moment 
presented  new  challenges.  When  fitting  distribution  parameters  to  nominal  hourly  values,  the 
fitted  distributions  initially  provided  poor  matches  to  computed  exceedence  probabilities.  These 
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poor  matches  occurred  because  SHIPM07  gives  non-zero  (but  very  small)  lateral  loads  in  head 
seas,  which  arise  from  a  non-zero  value  when  evaluating  the  cosine  of  180  degrees.  To  avoid 
the  resulting  problems  when  fitting  distributions  to  ranked  nominal  hourly  maximum  loads,  the 
fitting  procedures  exclude  nominal  hourly  maximum  loads  which  are  less  than  0.001  times  the 
highest  nominal  hourly  maximum. 


7.1  Vertical  Bending  Moment  at  Midships 


Figures  1  to  6  show  vertical  bending  moment  at  midships  versus  exceedence  probability.  For 
the  nominal  hourly  sea  loads,  an  hourly  exceedence  probability  of  a  =  0.10  was  used  for  the 
sample  case  as  it  was  found  to  produce  similar  maximum  lifetime  loads  for  the  nominal  hourly 
and  load  cycle  approaches.  The  nominal  seaway  approach  produces  significantly  smaller  lifetime 
loads  than  the  other  two  approaches.  These  smaller  lifetime  loads  are  due  to  the  smaller  average 
exposure  times  for  the  extreme  sea  conditions,  which  were  found  to  be  significantly  less  than 
one  hour.  In  reality,  an  extreme  sea  condition  will  either  occur  for  at  least  an  hour  or  will 
not  occur  during  the  ship  lifetime;  thus,  the  average  exposure  times  used  for  nominal  seaway 
maxima  produce  misleading  results.  Due  to  the  relatively  poor  results  of  the  nominal  seaway 
approach,  it  is  omitted  from  most  of  the  following  discussion. 

In  Figure  2,  the  plotted  lines  for  lifetime  vertical  bending  moment  from  the  nominal  hourly 
and  nominal  seaway  approaches  have  noticeable  discontinuities  because  the  probabilities  are 
computed  from  small  numbers  of  discrete  conditions.  The  fitted  distributions  to  nominal  hourly 
and  load  cycle  vertical  bending  moments  in  Figures  3  and  4  indicate  that  both  the  Gumbel  and 
Weibull  distributions  provide  good  fits  to  the  computed  data;  however,  the  Gumbel  distribu¬ 
tions  have  the  disadvantage  of  incorrectly  giving  negative  bending  moments  for  load  cycles  as 
exceedence  probability  approaches  1.0.  Consequently,  the  Weibull  distribution  is  recommended 
for  modelling  nominal  hourly  maximum  or  load  cycle  amplitude.  For  lifetime  extreme  bending 
moment,  Figures  5  and  6  indicate  that  the  Gumbel  distribution  provides  a  significantly  better 
fit  than  the  Weibull  distribution.  The  superior  fit  of  the  Gumbel  distribution  is  expected  for 
lifetime  maxima  because  it  is  based  on  the  asymptotic  behaviour  of  largest  extreme  values. 

Table  2  gives  means  and  standard  deviations  for  fitted  distributions  of  vertical  bending  mo¬ 
ment  at  midships.  The  nominal  hourly  and  load  cycle  amplitude  values  are  based  on  Weibull 
distributions,  while  the  lifetime  maxima  are  based  on  Gumbel  distributions.  It  should  be  em¬ 
phasized  that  the  properties  of  the  lifetime  maximum  based  on  nominal  hourly  maxima  will 
depend  on  the  hourly  exceedence  probability  a  used  to  determine  the  nominal  maximum.  For 
the  present  case  of  a  =  0.1,  the  lifetime  maxima  are  similar  for  the  two  different  approaches. 
The  lifetime  maximum  design  vertical  bending  moment  depends  on  the  probability  of  excee¬ 
dence  8  for  the  lifetime  maximum.  For  8  =  0.01,  the  design  vertical  bending  moment  would  be 
approximately  780  MN-m. 
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Figure  1:  Vertical  Bending  Moment  at  Midships  Versus  Nominal  Hourly,  Nominal  Seaway  and 
Load  Cycle  Exceedence  Probabilities 


Figure  2:  Vertical  Bending  Moment  at  Midships  Versus  Lifetime  Exceedence  Probability 
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Figure  3:  Computed  and  Fitted  Distributions  for  Nominal  Hourly  Vertical  Bending  Moment  at 
Midships,  a  —  0.1 


Figure  4:  Computed  and  Fitted  Distributions  for  Load  Cycle  Vertical  Bending  Moment  at 
Midships 
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Figure  5:  Computed  and  Fitted  Distributions  for  Lifetime  Vertical  Bending  Moment  at  Midships 
from  Nominal  Hourly  Loads,  a  =  0.1 


Figure  6:  Computed  and  Fitted  Distributions  for  Lifetime  Vertical  Bending  Moment  at  Midships 
from  Load  Cycles 
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Table  2:  Mean  and  Standard  Deviation  Values  for  Vertical  Bending  Moment  at  Midships 


Description 

Mean 

(MN-m) 

Standard  deviation 
(MN-m) 

Nominal  hourly,  a  =  0.1  (Weibull) 

66.1 

60.5 

Load  cycle  amplitude  (Weibull) 

22.5 

26.5 

Lifetime  maximum  from  hourly,  a  =  0.1  (Gumbel) 

641.2 

42.7 

Lifetime  maximum  from  cycle  (Gumbel) 

631.2 

47.2 

Figures  7  to  9  give  conditional  probabilities  for  ship  heading,  wave  period,  and  wave  height 
when  vertical  bending  moment  has  an  amplitude  greater  than  600  MN-m.  These  figures  indicate 
which  conditions  are  predominant  for  large  vertical  bending  moments.  Not  surprisingly,  the  large 
vertical  bending  moments  are  most  likely  to  be  caused  by  head  seas  and  large  waves.  The  most 
likely  wave  periods  for  large  vertical  bending  moments  are  associated  with  relatively  steep  waves. 


Figure  7:  Conditional  Probability  of  Heading  Given  Vertical  Bending  Moment  Amplitude 
Greater  than  600  MN-m 
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Significant  Wave  Height  Hs  (m) 


Figure  8:  Conditional  Probability  of  Significant  Wave  Height  Given  Vertical  Bending  Moment 
Amplitude  Greater  than  600  MN-m 


Zero-Crossing  Wave  Period  Tz  (s) 


Figure  9:  Conditional  Probability  of  Zero-Crossing  Wave  Period  Given  Vertical  Bending  Moment 
Amplitude  Greater  than  600  MN-m 
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7.2  Horizontal  Bending  Moment  at  Midships 

Figures  10  to  15  show  horizontal  bending  moment  at  midships  versus  exceedence  probability. 
The  general  trends  are  essentially  the  same  as  for  vertical  bending  moment.  The  lifetime  ex¬ 
ceedence  probabilities  are  essentially  the  same  for  the  nominal  hourly  maximum  and  load  cycle 
approaches  when  the  nominal  hourly  exceedence  probability  a  is  set  to  0.1.  As  was  observed  for 
vertical  bending  moment,  the  nominal  seaway  maximum  approach  predicts  significantly  smaller 
lifetime  loads  than  the  other  two  approaches. 

The  Weibull  distribution  provides  good  fits  to  exceedence  probabilities  for  nominal  hourly 
maximum  and  load  cycle  values;  however,  the  Gumbel  distribution  provides  superior  matching 
to  lifetime  exceedence  probabilities. 


Figure  10:  Horizontal  Bending  Moment  at  Midships  Versus  Nominal  Hourly  Maximum  and  Load 
Cycle  Exceedence  Probabilities 

Table  3  gives  means  and  standard  deviations  based  on  fitted  distributions.  The  values  of  the 
horizontal  bending  moment  parameters  are  approximately  one  third  of  the  values  for  vertical 
bending  moment.  Based  on  the  lifetime  maximum  parameters  and  a  probability  of  exceedence 
of  6  =  0.01,  the  design  horizontal  bending  moment  at  midships  would  be  approximately  235 
MN-m. 

Figures  16  to  18  give  conditional  probabilities  for  ship  heading,  wave  period,  and  wave  height 
when  horizontal  bending  moment  has  an  amplitude  greater  than  180  MN-m.  Largest  horizontal 
bending  moments  are  most  likely  to  occur  at  headings  of  60  and  120  degrees,  which  are  both 
30  degrees  off  beam  seas.  As  expected,  significant  wave  heights  are  likely  to  be  large  for  large 
bending  moments.  The  conditional  distribution  of  wave  periods  indicates  that  steep  waves  are 
dominant  for  high  horizontal  bending  moments. 
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Figure  11:  Horizontal  Bending  Moment  at  Midships  Versus  Lifetime  Exceedence  Probability 


Figure  12:  Computed  and  Fitted  Distributions  for  Nominal  Hourly  Maximum  Horizontal  Bend¬ 
ing  Moment  at  Midships,  a  =  0.1 
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Figure  13:  Computed  and  Fitted  Distributions  for  Load  Cycle  Horizontal  Bending  Moment  at 
Midships 


Figure  14:  Computed  and  Fitted  Distributions  for  Lifetime  Horizontal  Bending  Moment  at 
Midships  from  Nominal  Hourly  Loads,  a  =  0.1 
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Figure  15:  Computed  and  Fitted  Distributions  for  Lifetime  Horizontal  Bending  Moment  at 
Midships  from  Load  Cycles 


Table  3:  Mean  and  Standard  Deviation  Values  for  Horizontal  Bending  Moment  at  Midships 


Description 

Mean 

(MN-m) 

Standard  deviation 
(MN-m) 

Nominal  hourly  maximum,  a  =  0.1  (Weibull) 

17.3 

16.4 

Load  cycle  amplitude  (Weibull) 

6.0 

7.23 

Lifetime  maximum  from  hourly,  a  =  0.1  (Gumbel) 

183.4 

17.3 

Lifetime  maximum  from  cycle  (Gumbel) 

188.6 

14.1 
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Ship  Heading  0S  (degrees) 

Figure  16:  Conditional  Probability  of  Heading  Given  Horizontal  Bending  Moment  Amplitude 
Greater  than  180  MN-m 


Significant  Wave  Height  Hs  (m) 


Figure  17:  Conditional  Probability  of  Significant  Wave  Height  Given  Horizontal  Bending  Mo¬ 
ment  Amplitude  Greater  than  180  MN-m 
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Zero-Crossing  Wave  Period  Tz  (s) 


Figure  18:  Conditional  Probability  of  Zero-Crossing  Wave  Period  Given  Horizontal  Bending 
Moment  Amplitude  Greater  than  180  MN-m 
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8  Recommendations  for  Future  Work 


The  present  study  represents  an  initial  phase  for  computing  design  loads  for  ships,  and  can 
form  the  basis  for  more  sophisticated  analysis.  This  section  describes  several  areas  that  should 
be  considered  to  obtain  more  accurate  design  wave  loads. 

8.1  Wave  Scattergrams  for  Ship  Missions 

The  example  case  in  the  present  study  considers  a  ship  in  a  single  operational  area.  A  more 
realistic  analysis  should  use  a  composite  wave  scattergram  based  on  the  areas  and  seasons  in 
which  the  ship  will  be  operational.  Derivation  of  composite  scattergrams  is  a  relatively  simple 
procedure.  DREA  intends  to  develop  software  that  will  produce  composite  scattergrams  based 
on  BMT  Global  Wave  Statistics.  Saving  the  Global  Wave  Statistics  database  in  accessible  format 
is  the  biggest  task  of  this  development  and  is  currently  underway. 

8.2  Correlations  between  Wave  Conditions  and  Ship  Operations 

The  present  study  neglects  correlations  between  wave  conditions  and  ship  operations.  In 
reality,  ship  operations  can  be  greatly  influenced  by  wave  conditions,  particularly  in  severe  seas. 
For  example,  a  ship  captain  will  reduce  speed  in  severe  seas  to  reduce  the  incidence  of  slamming. 
Linear  wave  loads  appear  to  have  moderate  dependence  on  ship  speed  and  large  dependence  on 
ship  heading.  For  nonlinear  wave  loads  such  as  slamming,  they  can  be  greatly  dependent  on 
both  ship  speed  and  heading;  thus,  correlations  between  wave  conditions  and  ship  operations 
must  be  considered.  Ship  logs  can  be  used  to  determine  the  influence  of  wave  conditions  on  ship 
operations. 

8.3  Three  Dimensional  Hydrodynamic  Effects 

Load  computations  for  the  CPF  [7,  11]  have  shown  that  strip  theory  overpredicts  vertical 
bending  moment  at  midships  by  approximately  30  percent,  while  the  three  dimensional  code 
PRECAL  gives  better  results,  particularly  at  lower  ship  speeds.  The  wide  transom  stern  of  the 
CPF  necessitates  use  of  three  dimensional  hydrodynamic  coefficients  for  the  accurate  prediction 
of  sea  loads. 

8.4  Nonlinear  Wave  Loads 

Nonlinear  effects  can  significantly  influence  wave  loads,  particularly  for  severe  conditions 
which  are  of  greatest  interest  for  design.  Important  nonlinear  effects  include  non-zero  mean 
wave  loads,  sag/hog  ratios  greater  than  unity,  green  seas,  and  slamming.  Experiments  for  the 
CPF  reported  in  Reference  11  indicate  that  slamming  loads  are  very  important  but  that  non-zero 
mean  wave  loads  are  less  important. 

Guedes  Soares  and  Schellin  [3]  and  Jensen  and  Dogliani  [4]  give  frequency  domain  approaches 
for  predicting  wave  loads  based  on  nonlinear  strip  theory.  Alternatively,  Shin  et  al.  [12]  and 
Sclavounos  et  al.  [13]  give  time  domain  approaches  for  computing  nonlinear  sea  loads.  The  fre¬ 
quency  domain  approaches  have  the  advantages  of  being  computationally  efficient  and  directly 
producing  statistical  information.  If  a  nonlinear  frequency  domain  approach  were  to  be  used  for 
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the  CPF,  strip  theory  would  likely  be  inadequate  because  of  significant  three  dimensional  hydro- 
dynamic  effects.  Time  domain  approaches  are  likely  more  suitable  for  predicting  transient  loads 
from  Green  seas  and  slamming.  Due  to  difficulties  with  time  domain  predictions,  verification 
with  frequency  domain  predictions  and  experiments  is  usually  necessary. 

8.5  Wave  Data  for  Large  Wave  Heights 

The  sample  computations  in  the  previous  section  indicate  that  the  design  lifetime  sea  loads 
for  the  CPF  will  be  associated  with  significant  wave  heights  of  the  order  of  15  m.  Accurate 
probability  distributions  are  needed  for  wave  heights  associated  with  design  conditions;  how¬ 
ever,  BMT  Global  Wave  Statistics  has  the  limitation  of  a  single  wave  height  range  for  significant 
wave  heights  above  14  m.  Improved  modelling  of  wave  height  distributions  could  likely  be 
achieved  by  fitting  distributions  to  the  upper  wave  height  ranges  of  the  BMT  data.  Alterna¬ 
tively,  Reference  14  gives  distributions  for  wave  height  extremes  based  on  observed  storm  data. 
McTaggart  [15]  applied  extreme  wave  data  to  prediction  of  ship  capsize  in  waves. 

9  Conclusions 

Three  methods  have  been  successfully  developed  for  predicting  sea  load  exceedence  probabil¬ 
ities  based  upon  linear  theory.  Although  the  nominal  hourly  and  nominal  seaway  approaches  are 
useful,  the  load  cycle  approach  appears  to  be  superior  because  it  fully  considers  both  short-term 
and  long-term  variations  of  sea  loads. 

For  fatigue  computations,  the  load  cycle  approach  gives  required  exceedence  probabilities 
for  individual  load  cycles.  The  nominal  hourly,  nominal  seaway,  and  load  cycle  approaches  all 
give  lifetime  exceedence  probabilities  for  sea  loads;  however,  the  nominal  hourly  and  nominal 
seaway  load  approaches  require  that  the  user  select  an  exceedence  probability  for  hourly  and 
seaway  maximum  loads.  All  three  computational  approaches  can  yield  information  regarding 
conditions  associated  with  design  events. 

Sample  computations  for  the  CPF  demonstrate  that  the  approaches  can  be  used  to  determine 
probability  distributions  for  sea  loads.  The  nominal  seaway  approach  appears  to  significantly 
underpredict  lifetime  loads  because  the  expected  exposure  times  for  extreme  conditions  were 
found  to  be  short.  Gumbel  distributions  provide  very  good  fits  to  computed  exceedence  prob¬ 
abilities  for  lifetime  vertical  and  horizontal  bending  moments.  For  individual  load  cycles  and 
nominal  hourly  maxima,  Weibull  distribution  give  good  fits  to  computed  exceedence  probabil¬ 
ities.  For  maximum  lifetime  horizontal  and  vertical  bending  moments  for  the  CPF,  the  most 
likely  significant  wave  height  will  be  approximately  15  m. 

The  current  linear  analysis  could  form  the  basis  for  more  sophisticated  approaches.  Develop¬ 
ment  of  composite  wave  scattergrams  for  mission  and  lifetime  profiles  will  likely  be  an  easy  task. 
Correlations  between  wave  conditions  and  ship  operations  are  likely  very  important  for  severe 
conditions  and  should  be  considered.  Three  dimensional  hydrodynamic  effects  are  important 
for  the  CPF  and  can  be  included  using  codes  such  as  PRECAL.  Nonlinear  effects,  which  can 
be  very  important  in  severe  seas,  could  be  determined  using  nonlinear  codes  in  the  frequency 
or  time  domains.  While  time  domain  codes  offer  better  potential  for  transient  effects  such  as 
slamming  and  green  seas,  they  are  still  evolving  for  application  to  design. 


21 


A  EXTRMLIN  Input 

A.l  Input  Records 

Detailed  descriptions  of  EXTRMLIN  input  records  are  given  below.  Appendix  A. 2  gives 
sample  input  for  vertical  bending  moment  of  the  Canadian  Patrol  Frigate.  Each  new  input 
record  or  sub-record  corresponds  to  a  new  file  line.  The  format  of  the  input  file  may  be  adjusted 
by  inserting  extra  blanks  between  any  numerical  data,  and  by  placing  data  from  within  any 
particular  record  on  separate  lines;  however,  separate  records  cannot  be  combined  on  a  single 
line. 

Record  (a),  Eighty  Character  Title 
TITLE  (columns  1  -  80) 

TITLE  Alphanumeric  title  (maximum  of  80  characters)  which  is  written  on  output. 

Record  (b),  SHIPMOT  Post-processing  File 
SMPPRFILE  (columns  1-30) 

SMPPRFILE  Name  of  SHIPM07  binary  post-processing  file  (maximum  30  characters).  This 
file  must  have  metric  output  units  (OUTSYS  =  METRIC). 

Record  (c),  BMT  Global  Wave  Statistics  Wave  Scattergram  File 
WSFILE  (columns  1-30) 

WSFILE  Name  of  wave  scattergram  file  (maximum  30  characters).  The  scattergram  file 
has  joint  observations  of  significant  wave  height  and  wave  period  in  BMT  Global 
Wave  Statistics  ASCII  format. 

Record  (d),  Number  of  Ship  Speeds 

NSPEED  (1  integer) 

NSPEED  Number  of  ship  speeds.  This  must  correspond  to  the  number  of  ship  speeds 
NSPEED  in  the  SHIPM07  post-processing  file. 

Record  (dl),  Ship  Speeds  and  Probabilities 

SPDKNOT(I),  PSPEED(I)  (NSPEED  records  of  2  reals) 

SPDKNOT(I)  Ship  speed  in  knots.  This  must  correspond  to  ship  speed  SPDKNOT(I)  in  the 
SHIPM07  post-processing  file. 

PSPEED(I)  Probability  of  ship  travelling  at  speed  SPDKNOT(I).  The  PSPEED(I)  values 
should  have  a  sum  of  1 .0. 

Record  (e),  Number  of  Sea  Directions 

NSEADIR  (1  integer) 

NSEADIR  Number  of  sea  directions.  This  must  correspond  to  the  number  of  ship  speeds 
NSEADIR  in  the  SHIPM07  post-processing  file. 
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Record  (el),  Sea  Directions  and  Probabilities 

SEADIR(I),  PSEADIR(I)  (NSEADIR  records  of  2  reals) 

SEADIR(I)  Sea  direction  (from)  relative  to  ship  heading  (degrees).  This  must  correspond 
to  sea  direction  SEADIR(I)  in  the  SHIPM07  post-processing  file. 

PSEADIR(I)  Probability  of  waves  approaching  ship  from  direction  PSEADIR(I).  The  sum  of 
the  PSEADIR(I)  values  should  equal  1.0. 

Record  (f),  Exceedence  Probability  of  Nominal  Load  in  One  Hour  Seaway  and 

Variable  Duration  Seaway 
PMAXSEAWAY  (1  real) 

PMAXSEAWAY  Probability  of  exceedence  for  nominal  maximum  load  in  one  hour  seaway  and 
variable  duration  seaway. 


Record  (g),  Sea  Load  Mode 

LOADMODE  (1  character  string) 


LOADMODE  Control  string  for  sea  load  mode.  The  five  possible  inputs  are: 


HORSHEAR 

VERTSHEAR 

TORSION 

VERTBEND 

HORBEND 


Horizontal  shear  force. 
Vertical  shear  force. 
Torsional  moment. 

Vertical  bending  moment. 
Horizontal  bending  moment. 


Record  (h),  Station  for  Sea  Loads 

STATION  (1  real) 

STATION  Station  for  sea  loads.  The  station  must  correspond  to  a  station  XSTLOAD(I) 
in  the  SHIPM07  post-processing  file. 

Record  (i),  Ship  Life  and  Fraction  of  Life  Spent  at  Sea 
SHIPLIFE,  PATSEA  (2  reals) 

SHIPLIFE  Total  ship  life  (years). 

PATSEA  Fraction  of  ship  life  spent  at  sea. 

Record  (j),  Probability  Range  for  Fitting  Distributions 
PEXMINFIT,  PEXMAXFIT  (2  reals) 

PEXMINFIT  Lower  limit  on  exceedence  probability  for  fitting  Gumbel  and  Weibull  distribu¬ 
tions  to  computed  exceedence  probabilities. 

PEXMAXFIT  Upper  limit  on  exceedence  probability  for  fitting  Gumbel  and  Weibull  distribu¬ 
tions  to  computed  exceedence  probabilities. 
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Record  (k),  Increment  and  Maximum  of  Output  Load  Cycle  Amplitudes 

XCYCLEINC,  XCYCLEMAX  (2  reals) 

XCYCLEINC  Increment  of  output  load  cycle  amplitudes  (MN  or  MN-m). 

XCYCLEMAX  Maximum  output  load  cycle  amplitude  (MN  or  MN-m). 

Record  (1),  Lower  Threshold  of  Load  Cycle  Amplitude  for  Conditional  Probabilities 

XCY CLECOND  (1  real) 

XCYCLECOND  Lower  threshold  of  load  cycle  amplitude  for  conditional  probabilities  of  ship 
speed,  heading,  and  wave  conditions  (N  or  N-m). 


24 


A. 2  Sample  Input 

Midships  vertical  bending  moment  for  CPF 
cpfextreme.ppr 
areal 5 . gwa 
2 

10.0  0.5 
18.0  0.5 
7 

0.0  0.08325 
30.0  0.1667 
60.0  0.1667 
90.0  0.1667 
120.0  0.1667 
150.0  0.1667 
180.0  0.08325 
0.1 

VERTBEND 

10.0 

30.0  0.3 
0.0  0.9 
20.0  800.0 
600.0 


<  -  Record  (a)  Title 

<  -  Record  (b)  SHIPM07  file 

<  -  Record  (c)  Wave  statistics  file 

<  -  Record  (d)  Number  of  ship  speeds 

<  -  Record  (dl)  Speed,  probability 

<  -  Record  (dl) 

<  -  Record  (e)  Number  of  headings 

<  -  Record  (el)  Heading,  probability 

<  -  Record  (el) 

<  -  Record  (el) 

<  -  Record  (el) 

<  -  Record  (el) 

<  -  Record  (el) 

<  -  Record  (el) 

<  -  Record  (f)  Hourly  P(exceed) 

<  -  Record  (g)  Mode 

<  -  Record  (h)  Station  number 

<  -  Record  (i)  Life,  fraction  at  sea 

<  -  Record  (j)  Fit  probability  range 

<  -  Record  (k)  Cycle  inc  and  max 

<  -  Record  (1)  Conditional  threshold 


r 
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B  EXTRMLIN  Output 

B.l  Description  of  Output 

Appendix  B.2  gives  sample  computations  for  vertical  bending  moment  of  the  Canadian 
Patrol  Frigate.  The  output  begins  with  user  input  values,  which  are  followed  by  the  wave 
scattergram.  The  wave  scattergram  values  are  given  as  number  per  million  observations,  which 
corresponds  with  the  resolution  of  the  BMT  Global  Wave  Statistics  database.  Due  to  rounding 
error,  a  correction  factor  must  be  applied  to  the  wave  scattergram  data  to  ensure  that  the 
sum  of  probabilities  is  equal  to  one.  The  output  also  gives  information  from  the  SHIPM07 
post-processing  file. 

The  first  set  of  output  computations  is  for  nominal  hourly  maximum  loads.  The  output  gives 
nominal  hourly  sea  loads  and  conditions  for  ranked  events  associated  with  various  exceedence 
probabilities.  A  second  table  gives  exceedence  probabilities  and  conditions  for  the  highest  ranked 
load  events.  The  output  then  gives  statistical  fits  for  both  nominal  hourly  maximum  and  lifetime 
maximum  sea  loads.  After  the  nominal  hourly  sea  loads,  the  output  has  similar  computation 
results  for  the  nominal  seaway  sea  loads. 

The  third  set  of  output  computations  is  for  load  cycle  amplitudes.  Load  cycle  and  lifetime 
exceedence  probabilities  are  given  for  load  levels  requested  by  the  user.  Subsequent  output  gives 
statistical  fits  to  load  cycle  and  lifetime  maximum  loads.  The  output  finishes  with  conditional 
probabilities  for  ship  and  sea  conditions  when  loads  exceed  an  input  threshold  level. 
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B.2  Sample  Output 

Output  from  program  EXTRMLIN,  Risk  analysis  of  loads  using  linear  theory 
Defence  Research  Establishment  Atlantic 
Program  Version  1.0  -  January  1998 

EXTRMLIN  Run  Title: 

Midships  vertical  bending  moment  for  CPF 
11:27:23  08-May-98 


***************  ECHO  OF  USER  INPUT  *************** 

SHIPM07  post-processing  file  :  cpf extreme .ppr 

BMT  Global  Wave  Statistics  file  :  areal5.gwa 

2  ship  speeds 
Speed  (kt)  P (Speed) 

10.00  .500 

18.00  .500 


1.000 


7  ship  headings 

Heading  (deg) 

P (Heading) 

.00 

.083 

30.00 

.167 

60.00 

.167 

90.00 

.167 

120.00 

.167 

150.00 

.167 

180.00 

.083 

1.000 

Probability  of  exceedence  for  nominal  maximum  load  in  1  hour  seaway  : 
Load  mode  :  VERTBEND 

Ship  station  for  load  computation  :  10.0000 

Ship  life  (years)  :  30.00 

Fraction  of  life  at  sea  :  .3000 


Minimum  probability  of  exceedence  for  fitting  distributions  :  .000000E+00 
Maximum  probability  of  exceedence  for  fitting  distributions  :  .900000 

Load  cycle  parameters  (input  units  MN  or  MNm) 

Load  increment  :  20.000 

Maximum  load  :  800.000 

Conditional  design  load  :  600.000 


.1000 
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****************  Echo  of  BMT  Global  Wave  Statistics  Wave  Scattergram  **************** 


File  name  :  area!5.gwa 

Date  and  time  record  created  :  10/10/1997  8:39:29 


Sea  area  number 
Season  number 
Months 

Direction  number 
Percentage  of  obs  for  area 
Directions 
Direction  label 


15 

0 

January  to  December 
0 

100.000 
.0  to  360.0 
ALL  DIRECTIONS 


Number  of  observations  per  million 
I  Wave  Period  Tz  (s) 


Hs  (m)  | 

3.5 

4.5 

5.5 

6.5 

7.5 

8.5 

9.5 

10.5 

11.5 

12.5 

13.5 

.5  | 

195 

5148 

22478 

31077 

18074 

5512 

1059 

147 

16 

2 

1.5  I 

6 

1241 

20714 

71775 

86236 

49613 

16978 

4018 

730 

110 

14 

2.5  I 

230 

7349 

43695 

83798 

72651 

35743 

11681 

2838 

555 

93 

3.5  I 

45 

2203 

18702 

49097 

56350 

35688 

14653 

4377 

1034 

205 

4.5  I 

9 

631 

6984 

23184 

32851 

25176 

12298 

4304 

1175 

267 

5.5  I 

2 

181 

2480 

9923 

16639 

14855 

8342 

3316 

1018 

257 

6.5  I 

1 

54 

876 

4092 

7893 

8001 

5046 

2232 

756 

209 

7.5  | 

17 

316 

1683 

3659 

4136 

2881 

1396 

514 

154 

8.5  | 

5 

118 

704 

1700 

2113 

1605 

842 

334 

107 

9.5  | 

2 

46 

303 

803 

1085 

890 

502 

212 

72 

10.5  | 

1 

19 

135 

388 

566 

497 

299 

134 

48 

11.5  | 

8 

62 

192 

300 

281 

179 

85 

32 

12.5  1 

3 

29 

97 

162 

161 

108 

54 

21 

13.5  1 

2 

14 

51 

89 

93 

66 

34 

14 

14.5  | 

1 

15 

58 

115 

134 

105 

61 

28 

Correction  factor  applied  to  wave  scattergram  :  .999878 
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*******************  Echo  of  SHIPM07  Post -pro cessing  File  ******************* 


File  :  cpf extreme .ppr 

Title  :  Extreme  load  prediction  example  -  CPF 

16-APR-98  10:51:04 

Units  :  METRIC 
Spectrum  :  BRETSCHNEIDER 


Ship  characteristics 
Length 
Displacement 
LCB  aft  of  FP 
Draft  at  midships 
Trim  by  stern 
KG 

Beam  at  midships 


124.500  m 
4731.8  tonnes 
64.553  m 
5.025  m 
-.135  m 
6.230  m 
14.803  m 
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****  Nominal  hourly  maximum  loads  **** 


P (exceed)  Vertical  bending  moment  Ship  conditions  Wave  conditions 


Life 

Hour 

Max 

RMS 

Tz 

Speed 

Head 

Hs 

Tz 

(MNm) 

(MNm) 

(s) 

(kt) 

(deg) 

(m) 

(s) 

.057 

.000001 

703.1 

166.8 

5.0 

18.0 

150.0 

14.0 

6.5 

.563 

.000010 

619.9 

147.6 

5.3 

18.0 

180.0 

12.0 

7.5 

1.000 

.000103 

522.4 

124.9 

5.7 

18.0 

180.0 

12.0 

9.5 

1.000 

.000991 

390.2 

93.8 

6.2 

10.0 

150.0 

9.0 

7.5 

1.000 

.009721 

264.4 

67.1 

15.4 

10.0 

.0 

13.0 

13.5 

1.000 

.049942 

181.9 

43.6 

5.9 

10.0 

150.0 

4.0 

6.5 

1.000 

.100924 

145.3 

34.6 

5.3 

18.0 

150.0 

3.0 

7.5 

1.000 

.200888 

109.8 

28.5 

21.6 

18.0 

30.0 

4.0 

9.5 

1.000 

.299333 

88.4 

22.4 

15.0 

10.0 

.0 

3.0 

10.5 

1.000 

.400456 

71.8 

18.0 

12.5 

18.0 

60.0 

3.0 

9.5 

1.000 

.500437 

56.6 

13.5 

5.6 

10.0 

120.0 

2.0 

7.5 

1.000 

.597382 

41.7 

9.9 

5.1 

18.0 

90.0 

14.0 

9.5 

1.000 

.699032 

33.1 

7.9 

5.4 

10.0 

90.0 

12.0 

9.5 

1.000 

.800062 

12.7 

3.0 

4.7 

18.0 

90.0 

3.0 

7.5 

1.000 

.897142 

3.9 

.9 

5.7 

18.0 

90.0 

2.0 

12.5 

1.000 

.949297 

.0 

.0 

5.9 

10.0 

150.0 

.0 

6.5 

1.000 

.990254 

.0 

.0 

13.0 

10.0 

30.0 

.0 

7.5 
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****  Highest  nominal  hourly  maximum  loads  for  observed  conditions  **** 


P (exceed)  Vertical  bending  moment  Ship  conditions  Wave  conditions 


Life 

Hour 

Max 

RMS 

Tz 

Speed 

Head 

Hs 

Tz 

(MNm) 

(MNm) 

(s) 

(kt) 

(deg) 

(m) 

(s) 

.048 

.000001 

723.2 

172.2 

5.3 

18.0 

180.0 

14.0 

7.5 

.051 

.000001 

718.9 

170.6 

5.0 

18.0 

180.0 

14.0 

6.5 

.057 

.000001 

703.1 

166.8 

5.0 

18.0 

150.0 

14.0 

6.5 

.146 

.000002 

678.2 

161.4 

5.3 

18.0 

150.0 

14.0 

7.5 

.294 

.000004 

675.6 

161.2 

5.5 

18.0 

180.0 

14.0 

8.5 

.326 

.000005 

671.5 

159.9 

5.3 

18.0 

180.0 

13.0 

7.5 

.330 

.000005 

667.5 

158.4 

5.0 

18.0 

180.0 

13.0 

6.5 

.339 

.000005 

652.9 

154.9 

5.0 

18.0 

150.0 

13.0 

6.5 

.341 

.000005 

652.5 

156.5 

6.1 

10.0 

180.0 

14.0 

6.5 

.373 

.000006 

649.2 

156.2 

6.4 

10.0 

180.0 

14.0 

7.5 

.377 

.000006 

636.6 

152.5 

5.9 

10.0 

150.0 

14.0 

6.5 

.431 

.000007 

629.8 

149.9 

5.3 

18.0 

150.0 

13.0 

7.5 

.519 

.000009 

627.4 

149.7 

5.5 

18.0 

180.0 

13.0 

8.5 

.563 

.000010 

619.9 

147.6 

5.3 

18.0 

180.0 

12.0 

7.5 

.701 

.000015 

616.7 

147.1 

5.5 

18.0 

150.0 

14.0 

8.5 

.704 

.000015 

616.2 

146.2 

5.0 

18.0 

180.0 

12.0 

6.5 

.797 

.000020 

609.5 

145.7 

5.7 

18.0 

180.0 

14.0 

9.5 

.816 

.000021 

607.0 

145.8 

6.2 

10.0 

150.0 

14.0 

7.5 

.817 

.000022 

605.9 

145.3 

6.1 

10.0 

180.0 

13.0 

6.5 

.849 

.000024 

602.9 

145.4 

6.6 

10.0 

180.0 

14.0 

8.5 
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****  Fitted  distributions  of  nominal  hourly  maximum  load  **** 


Vertical  bending  moment 
Weibull  distribution 
Hourly  P(exceed)  =  exp(-(x/k)**beta) 
Weibull  k  (MNm)  :  68.2213 

Weibull  beta  :  1.0921 


Mean 

Std  deviation 


65.9863  (MNm) 
60.4876  (MNm) 


GUMBEL  distribution 

Hourly  P(exceed)  =  1  -  exp (-( exp ( -alpha (x  -  u))) 
GUMBEL  u  (MNm)  :  35.3468 

GUMBEL  alpha  (/MNm)  :  .0196 

Mean  :  64.7899  (MNm) 

Std  deviation  :  65.4231  (MNm) 


P (exceed)  Vertical  bending  moment 


Observed 

Weibull 

Gumbel 

(MNm) 

(MNm) 

(MNm) 

.000001 

695.8 

755.4 

740.1 

.000010 

622.8 

639.2 

622.6 

.000100 

525.0 

521.1 

505.2 

.001000 

390.2 

400.4 

387.7 

.010000 

264.4 

276.2 

270.0 

.050000 

181.8 

186.3 

186.9 

.100000 

145.4 

146.4 

150.1 

.200000 

109.9 

105.5 

111.9 

.300000 

88.4 

80.9 

87.9 

.400000 

71.8 

63.0 

69.6 

.500000 

56.6 

48.8 

54.0 

.600000 

41.5 

36.9 

39.8 

.700000 

33.1 

26.5 

25.9 

.800000 

12.7 

17.3 

11.1 

.900000 

3.9 

8.7 

-7.2 
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****  Fitted  distributions  of  lifetime  maximum  nominal  hourly  maximum  load  **** 


Vertical  bending  moment 
Weibull  distribution 
Hourly  P(exceed)  =  exp(-(x/k)**beta) 
Weibull  k  (MNm)  :  657.8682 

Weibull  beta  :  17.7704 


Mean 

Std  deviation 


638.4595  (MNm) 
44.3610  (MNm) 


GUMBEL  distribution 

Hourly  P(exceed)  =  1  -  exp(-(exp(-alpha(x  -  u))) 
GUMBEL  u  (MNm)  :  621.8691 

GUMBEL  alpha  (/MNm)  :  .0301 

Mean  :  641.0717  (MNm) 

Std  deviation  :  42.6683  (MNm) 


P (exceed)  Vertical  bending  moment 

Observed  Weibull  Gumbel 


(MNm) 

(MNm) 

(MNm) 

.050000 

720.4 

699.8 

720.7 

.100000 

688.3 

689.5 

696.7 

.200000 

677.1 

675.7 

671.8 

.300000 

674.8 

664.8 

656.2 

.400000 

633.6 

654.6 

644.2 

.500000 

627.9 

644.4 

634.1 

.600000 

619.0 

633.5 

624.8 

.700000 

616.8 

620.8 

615.7 

.800000 

609.1 

604.6 

606.0 

.900000 

580.2 

579.6 

594.1 

33 


****  Nominal  seaway  maximum  loads  **** 


P(exceed) 

Vertical 

bending  moment 

Ship  conditions 

Wave  conditions 

Exposure  time 

Life 

Seaway 

Max 

RMS 

Tz 

Speed 

Head 

Hs 

Tz 

(MNm) 

(MNm) 

(s) 

Ckt) 

(deg) 

(m) 

(s) 

(hours) 

.001 

.000001 

582.5 

161.4 

5.3 

18.0 

150.0 

14.0 

7.5 

.10 

.020 

.000010 

524.6 

135.3 

5.7 

18.0 

180.0 

13.0 

9.5 

.29 

.181 

.000103 

460.7 

122.7 

6.4 

10.0 

180.0 

11.0 

7.5 

.20 

.855 

.000991 

354.8 

103.3 

12.3 

18.0 

60.0 

10.0 

6.5 

.12 

1.000 

.009721 

261.4 

115.2 

27.0 

18.0 

.0 

12.0 

6.5 

.01 

1.000 

.049942 

189.6 

67.3 

22.4 

18.0 

30.0 

8.0 

5.5 

.03 

1.000 

. 100924 

157.1 

44.4 

12.7 

18.0 

60.0 

14.0 

13.5 

00 

1.000 

.200888 

119.4 

33.6 

6.0 

10.0 

120.0 

13.0 

13.5 

.09 

1.000 

.299333 

97.9 

20.9 

9.5 

10.0 

60.0 

6.0 

11.5 

14.67 

1.000 

.400456 

80.2 

17.1 

13.8 

10.0 

.0 

2.0 

5.5 

24.11 

1.000 

.500437 

62.1 

12.0 

12.5 

18.0 

60.0 

2.0 

9.5 

234.85 

1.000 

.597382 

43.9 

14.4 

4.6 

10.0 

90.0 

13.0 

6.5 

.01 

1.000 

.699032 

35.1 

7.1 

21.6 

18.0 

30.0 

1.0 

9.5 

111.55 

1.000 

.800062 

13.9 

2.6 

5.4 

10.0 

90.0 

4.0 

9.5 

165.42 

1.000 

.897142 

4.9 

1.0 

5.9 

10.0 

90.0 

2.0 

11.5 

18.65 

1.000 

. 949297 

.0 

.0 

5.9 

10.0 

150.0 

.0 

6.5 

204.19 

1.000 

. 990254 

.0 

.0 

13.0 

10.0 

30.0 

.0 

7.5 

118.76 

****  Highest  nominal  seaway  maximum  loads  for  observed  conditions  **** 


P(exceed) 

Vertical 

bending  moment 

Ship  conditions 

Wave  conditions 

Life 

Seaway 

Max 

RMS 

Tz 

Speed 

Head 

Hs 

Tz 

Exposure  time 

(MNm) 

(MNm) 

(s) 

(kt) 

(deg) 

(m) 

(s) 

Hours 

.001 

.000001 

608.4 

161.2 

5.5 

18.0 

180.0 

14.0 

8.5 

.19 

.001 

.000001 

586.9 

172.2 

5.3 

18.0 

180.0 

14.0 

7.5 

.05 

.001 

.000001 

582.5 

161.4 

5.3 

18.0 

150.0 

14.0 

7.5 

.10 

.004 

.000002 

581.9 

147.1 

5.5 

18.0 

150.0 

14.0 

8.5 

.38 

.009 

. 000004 

574.5 

145.7 

5.7 

18.0 

180.0 

14.0 

9.5 

.38 

.010 

. 000005 

559.9 

149.7 

5.5 

18.0 

180.0 

13.0 

8.5 

.17 

.010 

.000005 

541.8 

159.9 

5.3 

18.0 

180.0 

13.0 

7.5 

.05 

.010 

.000005 

541.6 

145.4 

6.6 

10.0 

180.0 

14.0 

8.5 

.19 

.010 

.000005 

540.0 

138.2 

5.5 

18.0 

180.0 

12.0 

8.5 

.32 

.011 

.000006 

538.0 

149.9 

5.3 

18.0 

150.0 

13.0 

7.5 

.09 

.012 

.000006 

537.5 

130.5 

5.6 

18.0 

150.0 

14.0 

9.5 

.76 

.014 

.000007 

535.9 

136.6 

5.5 

18.0 

150.0 

13.0 

8.5 

.34 

.018 

.000009 

530.8 

147.6 

5.3 

18.0 

180.0 

12.0 

7.5 

.10 

.020 

.000010 

524.6 

135.3 

5.7 

18.0 

180.0 

13.0 

9.5 

.29 

.029 

.000015 

524.0 

156.2 

6.4 

10.0 

180.0 

14.0 

7.5 

.05 

.030 

.000015 

523.9 

138.3 

5.3 

18.0 

150.0 

12.0 

7.5 

.19 

.039 

. 000020 

519.5 

145.8 

6.2 

10.0 

150.0 

14.0 

7.5 

.10 

.041 

.000021 

517.0 

132.1 

6.5 

10.0 

150.0 

14.0 

8.5 

.38 

.041 

. 000022 

516.7 

126.7 

5.5 

18.0 

180.0 

11.0 

8.5 

.63 

.046 

.000024 

515.0 

129.6 

5.9 

18.0 

180.0 

14.0 

10.5 

.44 
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****  Fitted  distributions  of  nominal  seaway  maximum  load  **** 


Vertical  bending  moment 
Weibull  distribution 
Hourly  P(exceed)  =  exp(-(x/k)**beta) 
Weibull  k  (MNm)  :  71.7517 

Weibull  beta  :  1.1460 


Mean 

Std  deviation 


68.3556  (MNm) 
59.7931  (MNm) 


GUMBEL  distribution 

Hourly  P(exceed)  =  1  -  exp (~ (exp ( -alpha (x  -  u))) 
GUMBEL  u  (MNm)  :  48.6501 

GUMBEL  alpha  (/MNm)  :  .0223 

Mean  :  74.5894  (MNm) 

Std  deviation  :  57.6376  (MNm) 


P (exceed)  Vertical  bending  moment 


Observed 

Weibull 

Gumbel 

(MNm) 

(MNm) 

(MNm) 

.000001 

582.3 

709.4 

669.5 

.000010 

527.0 

605.1 

566.0 

.000100 

460.7 

498.0 

462.6 

.001000 

354.8 

387.5 

359.1 

.010000 

261.4 

272.0 

255.4 

.050000 

189.5 

186.9 

182.1 

.100000 

157.1 

148.6 

149.8 

.200000 

119.9 

108.7 

116.1 

.300000 

97.9 

84.4 

95.0 

.400000 

80.3 

66.5 

78.8 

.500000 

62.1 

52.1 

65.1 

.600000 

43.8 

39.9 

52.6 

.700000 

35.0 

29.2 

40.3 

.800000 

13.9 

19.4 

27.3 

.900000 

4.7 

10.1 

11.2 

35 


****  Fitted  distributions  of  lifetime  maximum  nominal  seaway  maximum  load  **** 


Vertical  bending  moment 
Weibull  distribution 
Hourly  P(exceed)  =  exp(-(x/k)**beta) 
Weibull  k  (MNm)  :  431.9839 

Weibull  beta  :  8.0046 


Mean 

Std  deviation 


406.8290  (MNm) 
60.3287  (MNm) 


GUMBEL  distribution 

Hourly  P(exceed)  =  1  -  exp (-(exp (-alpha (x  -  u))) 
GUMBEL  u  (MNm)  :  388.2857 

GUMBEL  alpha  (/MNm)  :  .0272 

Mean  :  409.4846  (MNm) 

Std  deviation  :  47.1045  (MNm) 


P (exceed)  Vertical  bending  moment 


Observed 

Weibull 

Gumbi 

(MNm) 

(MNm) 

(MNm) 

010000 

541.7 

522.8 

557.2 

050000 

507.3 

495.4 

497.4 

100000 

474.2 

479.4 

470.9 

200000 

456.7 

458.4 

443.4 

300000 

432.1 

442.1 

426.1 

400000 

417.2 

427.3 

413.0 

500000 

397.3 

412.7 

401.7 

600000 

392.4 

397.2 

391.5 

700000 

381.0 

379.8 

381.5 

800000 

362.7 

358.2 

370.8 

900000 

352.8 

326.1 

357.7 
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*****  Fatigue  load  cycle  computations  **** 


Vertical  bending  moment 

Period  for  mean  zero-crossing  frequency:  7.18  seconds 

Number  of  cycles  during  ship  life  :  .395E+08 


Load 

P(exceed) 

Cycles  during 

(MNm) 

Cycle 

Life 

20.000 

.425933 

1.000000 

. 168E+08 

40.000 

.186885 

1.000000 

. 738E+07 

60.000 

.083012 

1.000000 

. 328E+07 

80.000 

.038180 

1.000000 

. 151E+07 

100.000 

.018266 

1.000000 

. 722E+06 

120.000 

.009067 

1.000000 

. 358E+06 

140.000 

.004652 

1.000000 

. 184E+06 

160.000 

.002460 

1.000000 

. 972E+05 

180.000 

.001336 

1.000000 

. 528E+05 

200.000 

.000743 

1.000000 

. 293E+05 

220.000 

.000422 

1.000000 

. 167E+05 

240.000 

.000244 

1.000000 

. 963E+04 

260 . 000 

.000143 

1.000000 

. 565E+04 

280.000 

.000085 

1.000000 

. 336E+04 

300.000 

.000051 

1.000000 

. 202E+04 

320.000 

.000031 

1.000000 

. 122E+04 

340.000 

.000019 

1.000000 

. 747E+03 

360.000 

.000012 

1.000000 

.458E+03 

380.000 

. 71E-05 

1.000000 

. 282E+03 

400.000 

. 44E-05 

1.000000 

. 174E+03 

420.000 

. 27E-05 

1.000000 

. 107E+03 

440.000 

. 17E-05 

1.000000 

. 662E+02 

460.000 

. 10E-05 

1.000000 

.408E+02 

480.000 

. 64E-06 

1.000000 

.251E+02 

500.000 

. 39E-06 

1.000000 

. 154E+02 

520.000 

. 24E-06 

.999920 

. 943E+01 

540.000 

. 15E-06 

.996816 

. 575E+01 

560.000 

. 88E-07 

.969453 

. 349E+01 

580 . 000 

. 53E-07 

. 878302 

. 211E+01 

600.000 

. 32E-07 

.717722 

. 126E+01 

620.000 

.19E-07 

.530104 

. 755E+00 

640.000 

. 11E-07 

.361222 

. 448E+00 

660.000 

. 67E-08 

.232232 

. 264E+00 

680.000 

. 39E-08 

.143383 

. 155E+00 

700.000 

. 23E-08 

.086066 

. 900E-01 

720.000 

. 13E-08 

.050624 

. 520E-01 

740.000 

.75E-09 

.029323 

. 298E-01 

760.000 

.43E-09 

.016775 

. 169E-01 

780.000 

. 24E-09 

.009493 

. 954E-02 

800.000 

. 14E-09 

.005320 

. 533E-02 
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****  Fitted  distributions  of  load  cycle  amplitude 
Vertical  bending  moment 


Weibull  distribution 


P(exceed)  =  exp(-(x/k)**beta) 


Weibull  k  (MNm) 
Weibull  beta 
Mean 

Std  deviation 


20.8062 

.8536 

22.5766  (MNm) 
26.5549  (MNm) 


Gumbel  distribution 

P(exceed)  =  1  -  exp (- (exp ( -alpha (x  -  u))) 
Gumbel  u  (MNm)  :  -53.0322 


Gumbel  alpha  (/MNm)  : 

.0267 

Mean 

:  -31. 

4332  (MNm) 

Std  deviation 

:  47.9934  (MNm) 

Load 

P (Exceed) 

(MNm) 

Observed 

Weibull 

Gumbel 

20.00 

.425933 

.380288 

.132411 

40.00 

.186885 

.174281 

.079861 

60.00 

.083012 

.084618 

.047601 

80.00 

.038180 

.042552 

.028175 

100.00 

.018266 

.021939 

.016607 

120.00 

.009067 

.011532 

.009765 

140.00 

.004652 

.006156 

.005734 

160.00 

.002460 

.003330 

.003364 

180.00 

.001336 

.001821 

.001973 

200.00 

.000743 

.001006 

.001156 

220.00 

.000422 

.000560 

.000678 

240.00 

.000244 

.000315 

.000397 

260.00 

.000143 

.000178 

.000233 

280.00 

.000085 

.000101 

.000136 

300 . 00 

.000051 

.000058 

.000080 

•  320.00 

.000031 

.000033 

.000047 

340.00 

.000019 

.000019 

.000027 

360.00 

.000012 

.000011 

.000016 

380.00 

. 71E-05 

.65E-05 

. 94E-05 

400.00 

.44E-05 

. 38E-05 

. 55E-05 

420.00 

. 27E-05 

. 23E-05 

. 32E-05 

440.00 

. 17E-05 

. 13E-05 

. 19E-05 

460.00 

. 10E-05 

.79E-06 

. 11E-05 

480.00 

. 64E-06 

.47E-06 

. 65E-06 

500.00 

.39E-06 

. 28E-06 

. 38E-06 

520.00 

.24E-06 

. 17E-06 

. 22E-06 

540.00 

. 15E-06 

. 10E-06 

. 13E-06 

560.00 

.88E-07 

. 60E-07 

. 77E-07 

580.00 

.53E-07 

. 36E-07 

.45E-07 

600.00 

. 32E-07 

. 22E-07 

. 26E-07 

**** 
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620.00 

640.00 

660.00 

680.00 

700.00 

720.00 

740.00 

760.00 

780.00 

800.00 


. 19E-07 
.11E-07 
.67E-08 
.39E-08 
.23E-08 
. 13E-08 
. 75E-09 
.43E-09 
.  24E-09 
.  14E-09 


.  13E-07 
.  81E-08 
.49E-08 
.30E-08 
.18E-08 
.  11E-08 
.70E-09 
.43E-09 
.26E-09 
. 16E-09 


.  15E-07 
. 91E-08 
.53E-08 
.31E-08 
. 18E-08 
.11E-08 
.63E-09 
.  37E-09 
.  21E-09 
. 13E-09 
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****  Fitted  distributions  of  maximum  lifetime  load  cycle  amplitude  **** 


Vertical  bending  moment 


Weibull  distribution 


P(exceed)  =  exp(-(x/k)**beta) 


Weibull  k  (MNm) 
Weibull  beta 
Mean 

Std  deviation 


660.3448 

11.3893 

631.5430  (MNm) 
67.1749  (MNm) 


Gumbel  distribution 

P(exceed)  =  1  -  exp(-(exp(-alpha(x  -  u))) 
Gumbel  u  (MNm)  :  609.7956 


Gumbel  alpha  (/MNm)  : 

.0272 

Mean 

:  631. 

0432  (MNm) 

Std  deviation 

:  47. 

2125  (MNm) 

Load 

P(Exceed) 

(MNm) 

Observed  Weibull 

Gumbel 

20.00 

1.000000 

1.000000 

1.000000 

40.00 

1.000000 

1.000000 

1.000000 

60.00 

1.000000 

1.000000 

1.000000 

80.00 

1.000000 

1.000000 

1.000000 

100.00 

1.000000 

1.000000 

1.000000 

120.00 

1.000000 

1.000000 

1.000000 

140.00 

1.000000 

1.000000 

1.000000 

160.00 

1.000000 

1.000000 

1.000000 

180.00 

1.000000 

1.000000 

1.000000 

200.00 

1.000000 

.999999 

1.000000 

220.00 

1.000000 

.999996 

1.000000 

240.00 

1.000000 

.999990 

1.000000 

260.00 

1.000000 

.999975 

1.000000 

280.00 

1.000000 

.999943 

1.000000 

300.00 

1.000000 

.999875 

1.000000 

320.00 

1.000000 

.999739 

1.000000 

340.00 

1.000000 

.999479 

1.000000 

360.00 

1.000000 

.999002 

1.000000 

380.00 

1.000000 

.998154 

1.000000 

400.00 

1.000000 

.996691 

1.000000 

420.00 

1.000000 

.994239 

1.000000 

440 . 00 

1.000000 

. 990234 

1.000000 

460.00 

1.000000 

.983849 

1.000000 

480.00 

1.000000 

. 973907 

1.000000 

500.00 

1.000000 

.958784 

1.000000 

520.00 

.999920 

.936327 

.999990 

540.00 

.996816 

.903824 

.998718 

560.00 

.969453 

.858119 

.979099 

580.00 

.878302 

.795973 

.894240 

600.00 

.717722 

.714819 

.728792 
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620.00 
640.00 
660.00 
680.00 
700.00 
720.00 
740.00 
760.00 
780 . 00 
800.00 


.530104 

.361222 

.232232 

.143383 

.086066 

.050624 

.029323 

.016775 

.009493 

.005320 


.614023 
.496495 
.370068 
.247432 
. 143282 
.068705 
. 025767 
.007033 
.001276 
.000138 


.531349 

.356095 

.225611 

.138004 

.082640 

.048865 

.028679 

.016759 

.009769 

.005685 
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****  Conditional  Probabilities  for  Loads  Exceeding  Threshold  **** 

Conditional  probabilities  are  for  cycle  loads  exceeding  threshold  level 
Vertical  bending  moment 

Threshold  level  :  600.000  (MNm) 

Cycle  exceedence  probability  :  .320199E-07 


Condition  probabilities  for  ship  speed 
Ship  speed  P (speed) 

(knots) 


10.0  .1463 

18.0  .8537 


1.0000 


Condition  probabilities  for  ship  heading 
Ship  heading  P (heading) 

(deg) 


.0 

30.0 

60.0 

90.0 

120.0 

150.0 

180.0 


.0062 

.0041 

.0012 

.0000 

.0001 

.3565 

.6318 


1.0000 


Condition  probabilities  for  significant  wave  height 
Wave  height  P(Hs) 


.5 

.0000 

1.5 

.0000 

2.5 

.0000 

3.5 

.0000 

4.5 

.0000 

5.5 

.0000 

6.5 

.0000 

7.5 

.0000 

8.5 

.0002 

9.5 

.0021 

10.5 

.0130 

11.5 

.0451 

12.5 

.1057 

13.5 

.1994 

14.5 

.6345 

1.0000 
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Condition  probabilities  for  zero-crossing  wave  period 
Wave  period  P(Tz) 


3.5 

.0000 

4.5 

.0000 

5.5 

.0000 

6.5 

.0478 

7.5 

.3808 

8.5 

.4257 

9.5 

.1314 

10.5 

.0138 

11.5 

.0005 

12.5 

.0000 

13.5 

.0000 

1.0000 

Conditional  probabilities  for  wave  height  and  period  combinations 
Cell  values  are  number  of  occurences  per  million 
Blank  cells  indicate  zero  probability 


i 

Wave  Period  Tz 

(s) 

Hs  (m) 

i 

3.5  4.5 

5.5 

6.5 

.5 

i 

1.5 

i 

2.5 

i 

3.5 

i 

4.5 

i 

5.5 

i 

6.5 

i 

7.5 

i 

8.5 

i 

23 

9.5 

i 

3 

278 

10.5 

i 

20 

1449 

11.5 

i 

4159 

12.5 

i 

6938 

13.5 

i 

15142 

14.5 
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